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Enhanced stimulated Raman scattering in slow-light photonic crystal waveguides 

James F. McMillan, Xiaodong Yang, Nicolae C. Panoiu, Richard M. Osgood, and Chee Wei Wong^ 
'Optical Nano structures Laboratory, Columbia University, New York, NY 10027 

We investigate for the first time the enhancement of the stimulated Raman scattering in slow-light 
Silicon-on-Insulator (SOI) photonic crystal line defect waveguides. By applying the Bloch-Floquet 
formalism to the guided modes in a planar photonic crystal, we develop a formalism that relates the 
intensity of the down-shifted Stokes signal to the pump intensity and the modal group velocities. The 
formalism is then applied to two prospective schemes for enhanced stimulated Raman generation in 
slow-light photonic crystal waveguides. The results demonstrate a maximum factor of 10* (66,000) 
enhancement with respect to SOI channel waveguides. Effects of two photon absorption, intrinsic 
scattering, and disorder with respect to slow-light Raman generation towards optically-pumped 
silicon amplifiers and lasers are also discussed. 



Silicon photonics has seen remarkable advancements in 
recent years. Subwavelength silicon nanostructures - such 
as photonic crystals and high-index-contrast photonic in- 
tegrated circuits - offer the opportunity to manipulate 
the propagation of light at sub- wavelength scales. More- 
over, the inherent ease of integrating the silicon photonics 
platform with CMOS foundry ICs offers unprecedented 
bandwidth per unitcost and distance in optical data com- 
munications. 

Silicon, however, is at an intrinsic disadvantage for op- 
tical amplification and lasing due to its indirect band 
gap and non-existent second order nonlinear response. 
Recent work has demonstrated that stimulated Raman 
scattering (SRS) in single-crystal silicon channel waveg- 
uides is a feasible means to achieve amplification and 
lasing via optical pumpingiiSiSi^. This is due to the 
intrinsically large Raman gain coefRcicnt in silicon (be- 
ing 10"^ to 10'' times greater than for silica) and sili- 
con nanostructures offering the benefit of high optical 
confinement due to high-index contrast of silicon with 
air or silicon oxide. While still requiring an optical 
pump and possessing limited gain bandwidth, enhanced 
SRS through slow-light silicon photonic crystal waveg- 
uides (PhCWG) can serve as an ultra-compact on-chip 
gain media at desired telecommunications frequencies. 
Enhanced Raman scattering has been observed in bulk 
hollow-core slow-light guided-wave structures^ and has 
also recently been suggested for photonic crystal (PhC) 
defect nanocavitiesii In addition, a semiclassical model 
of Raman scattering in bulk photonic crystals has been 
introducedi^ In this Letter we demonstrate theoretically 
for the first time the explicit enhancement of SRS in a 
slow-light PhCWG through a four wave mixing formal- 
ism from the computed modes of the line-defect waveg- 
uide. A silicon PhCWG studied here, made by removing 
a single row in a hexagonal lattice of holes - denoted as 
"Wl PhCWG" - and its projected band structure can be 
seen in Fig. This structure supports two tightly con- 
fined modes with small group velocities, as illustrated 
by the two bands within the band gap, with frequen- 
cies below the light line. The field distribution of these 
two modes, as computed through the plane wave expan- 



sion method,^ is illustrated in Fig. [3 The strong sub- 
wavelength modal confinement of the high index contrast 
PhCWG leads to increased field intensities in the silicon 
gain media, permitting increased nonlinear interactions. 
In addition to increased field intensities from high in- 
dex confinement, there is additional SRS enhancement 
from the small group velocities of the PhCWG propagat- 
ing modes. Physically this enhancement originates from 
the effective long light-matter interaction times at small 
group velocities. Photon localization is observed at the 
band edge; the photon experiences multiple scattering 
processes and moves very slowly through the material 
structure. The guided bands of a 2D PhCWG can be de- 
signed to be as flat as desired {vg = div/dk) for slow-light 
behavior, and group velocities as low as 10~^c to 10~'^c 
have been demonstrated>i2iii 

In SRS for silicon an incident photon interacts with the 
LO and TO phonons. The strongest Stokes peak arises 
from the single first-order Raman-phonon at the center of 
the Brillouin zone. The generation of the Stokes photons 
can be understood classically as a third order nonlinear 
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FIG. 1: Projected band structure of silicon Wl PhCWG 
indicating pump and Stokes frequencies, (top) Scheme 1, 
(r/a = 0.29) (bottom) Scheme 2, {r/a = 0.22). {h/a = 0.6) 
in both cases, (inset) Wl PhCWG. 
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effect and this formalism has been used to model SRS 
in SOI waveguides, both in CW^^ and pulsedi^^ opera- 
tion. It can be modeled in bulk materials as a degen- 
erate four-wave-mixing problem involving the pump and 
Stokes beams. The important material parameter is the 
third order nonlinear Raman susceptibility, ■ For sili- 
con, at resonance, is defined by the components X^jij 
= -ix^ = -ill.2 X 10-18 m^ • {i,j = 1,2,3). An ad- 
ditional symmetry, imposed by the crystal point group 
{m3m for Si), is Xujj = O-^xf^y- These components, 
and their permutations as defined by the crystal point 
group, define the SRS in a silicon crystal. For our pur- 
pose we shall consider scattering in silicon along the [110] 
direction since practical devices are fabricated along this 
direction due to the favorable cleaving of silicon along 
this direction. 

For bulk silicon, the evolution of the Stokes beam is 
defined by the following equationiS 



dz 



3^Jm(xfff) 



(1) 



where XcS = Y^tjuXfju^Wihai- Here, a and P are 
unit vectors along the polarization directions of the pump 
and Stokes beams, respectively. Eq. Q describes the 
gain of the Stokes intensity. Is- It shows an intrinsic 
dependence on the polarization and the phonon selection 
rules through and the intensity of the pump beam 
by Ip. The bulk solution also describes SRS in dielectric 
waveguides, where x^s is averaged over the waveguide 
mode field distribution. 

A PhCWG presents a very different field distribution 
than the bulk or dielectric waveguide case. As shown in 
the computed modal profiles of Fig. the mode differs 
from that of a conventional channel waveguide in that it 
exhibits a periodic variation in the direction of propaga- 
tion. We introduce the modal distribution of the pump 
and Stokes modes in a Bloch-Floquct formalism, 

£nM^ (r, ^n) = E„,k„ (r, ^n) exp[ik(w„) • r], (2) 

where n is a mode index (n = p, s), k„ — k([j„) is the 
mode wave- vector, E„^k„(r,w„) is the modal distribu- 
tion within a unit cell of the PhC, defined in Fig. 12 
and obeys Bloch boundary condition E„.k„ (r + A, aj„) — 
E„,k„(r,w„). A defines the length of the unit cell in 
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FIG. 2: Calculated bound states of a hexagonal lattice Wl 
PhCWG with defect modes separated by the LO/TO optical 
phonon (Scheme 1). a)Stokes b)Pump 



the direction of propagation and for a Wl waveguide 
this equals the PhC lattice constant a. To develop an 
equation that relates the evolution of the Stokes mode 
to the pump mode, we employ the Lorentz reciprocity 
theorem, 1^ 

d 



>Z J A 



I A 

(3) 

This relates the unperturbed PhCWG modes of the 
pump or Stokes wavelengths, {E„k„,H„k„}: to those 
of the nonlinearly induced fields. The envelopes of the 
fields are defined as 

E(r) = Us(z)Es,k,,(r,t^s) + Mp(2;)Ep,kp(r,t^p), (4a) 
H(r) ^ Us{z)lisMA^^^s) + Up{z)'Rp,i,^{r,ujp), (4b) 

with the assumption that the change in the pump and 
Stokes field amplitudes, Up{z) and Us{z) respectively, 
over the length of the unit cell of the waveguide is very 
small (A^^^j^ <C 1). Taking the fields as defined in Eq. 
llHl, we derive the dependence of the Stokes amplitude on 
the longitudinal distance, z, 



dz 
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/ P«(r,c^,)-E,.k,(r,c^,)dF, (5) 



where Pg is the mode power and P^^ir, LUs) = 

6eoX^:E*_kX'')Ep.kXr)Es,k,(r)|uppUs. The integral in 
Eq. (O is taken over the volume (Vq) of the unit cell of 
the PhCWG mode. Furthermore, the group velocity of 
the modes can expressed by the following equationi^ 
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With Eqs. and and by rewriting Eq. (jSJ in terms 
of the modes intensity, an equation for the intensity of 
the Stokes mode inside the PhCWG is obtained. 



dh 
dz 



where 



rj I 

eovPgVg 



(7) 



A^efflni(^ E*(u;,) • x'':E*(wp)E(u;p)E(c^,)dI/ 

i 1^ e(r)|E(u;p)|W)Qj|^ e(r)|E(c..)|W 

(8) 

is the efffcctive susceptibility. Here, the effective area A^f 
is defined as the average modal area across the volume 
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The final equation, Eq. 0, shows the explicit inverse 
dependence the Stokes mode amplification has on the 
group velocities of the pump and Stokes modes. 

Table ^ shows the results of Eq. iQ as being applied 
to two different PhCWG schemes for SRS. The group 
velocities are calculated from the slope of the projected 
band structure. The first (Scheme 1) involves utilizing 
both the guided modes of the Wl waveguide; odd-parity 
is the pump mode and even-parity is the Stokes mode. 
The wavelength separation of the modes at the edge of 
the Brillioun zone is matched to the LO/TO frequency 
separation of the pump and Stokes beams (15.6 THz in 
Sii&). The second (Scheme 2) utilizes a wide bandwidth 
PhCWGfii in order for the Stokes and pump modes to 
exist both in the fundamental mode and below the light 
line. The arrows in Fig. |21 indicate the pump and Stokes 
frequency locations for both schemes. 

From the results of Tabled the Raman gain - propor- 
tional to K/v^Vg - is enhanced by up to approximately 
10^ (Scheme 1:66,000, Scheme 2:86) times compared to 
bulk Si based on a comparison of the respective group 
velocities. The results in Table |2 also show a k value 
of the same order with a conventional SOI waveguide^ 
In addition, we note a reduction in k in Scheme 1 as 
compared to Scheme 2, due to the lower modal overlap. 
However, the single mode (Scheme 2) operation has the 
disadvantage that only the Stokes mode, and not both 
modes, are at low group velocities for enhanced SRS. 

The above results highlight the benefits of SRS en- 
hancement through slow-light interactions in compact 
PhCWG schemes. This approach can be readily ex- 
tended to include two photon and bulk free carrier ab- 
sorption effects^'* which may limit the effective Raman 
gain in PhCWGs. These effects, in the experimental re- 
alization of silicon SRS amplification and lasing in slow- 
light PhCWGs, can be surmounted with pulsed-laser 
operation** or PIN diodes^ to sweep the free-carriers. 



In addition, we note recent theoreticalia and 
experimental^^ studies of PhCWGs, which show 
that slow group velocity modes exhibit increased 
scattering losses. These losses are from coupling and in- 
trinsic (backscatter) reflection. Coupling into slow-light 
modes is currently the dominant loss experimentally, 
although this can in principle be reduced through 
careful adiabatic coupling-^" between the PhCWGs and 
input /output channel bus waveguides. Moreover, with 
thorough attention to fabrication disorder, reflection 
losses in PhCWG are suggested to be comparable with 
index- guided waveguide^. These scattering losses can 
thus potentially be smaller than the enhanced SRS 



TABLE I: Group velocity and Effective Susceptibility 
in PhCWG schemes 



Scheme 


^« 


"9 


Ac(xlO-i^)[m2.V-2] 


1 


0.00017c 


0.0077c 


0.55 


2 


0.0041c 


0.24c 


2.02 



gain discussed, permitting the possibility for compact 
silicon Raman ampliflers and lasers. We also note that, 
for the same desired Raman gain, the device length 
is significantly reduced, by (c/wg)^, allowing compact 
integration for high-density photonic circuits. 
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